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ABSTRACT 

This paper reports a blind search for magnetospheric emissions from planets around nearby 
stars. Young stars are likely to have much stronger stellar winds than the Sun, and because 
planetary magnetospheric emissions are powered by stellar winds, stronger stellar winds may 
enhance the radio luminosity of any orbiting planets. Using various stellar catalogs, we selected 
nearby stars (< 30 pc) with relatively young age estimates (< 3 Gyr). We constructed different 
samples from the stellar catalogs, finding between 100 and several hundred stars. We stacked 
images from the 74- MHz (4-m wavelength) VLA Low- frequency Sky Survey (VLSS), obtaining 
3tT limits on planetary emission in the stacked images of between 10 and 33 mJy. These flux 
density limits correspond to average planetary luminosities less than 5-10 xlO^^ erg s~^. Using 
recent models for the scaling of stellar wind velocity, density, and magnetic field with stellar 
age, we estimate scaling factors for the strength of stellar winds, relative to the Sun, in our 
samples. The typical kinetic energy carried by the stellar winds in our samples is 15-50 times 
larger than that of the Sun, and the typical magnetic energy is 5-10 times larger. If we assume 
that every star is orbited by a Jupiter-like planet with a luminosity larger than that of the Jovian 
decametric radiation by the above factors, our limits on planetary luminosities from the stacking 
analysis are likely to be a factor of 10-100 above what would be required to detect the planets in 
a statistical sense. Similar statistical analyses with observations by future instruments, such as 
the Low Frequency Array (LOFAR) and the Long Wavelength Array (LWA), offer the promise 
of improvements by factors of 10-100. 

Subject headings: planetary systems 
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1. Introduction 

In searches for extrasolar planets, the precision 
to which a planetary signal can be extracted from 
the data can depend in part on the properties 
of the host star. For instance, in radial velocity 
surveys, one of the limiting factors in the veloc- 
ity precision is intrinsic stellar "jitter," caused 
by starspots or other surface inhomogeneities. 
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Such stellar jitter is well known to be corre- 
lated with st ellar activity, the lev el of which de- 
clines in age ( Butler et al. il996; Saar fc Donahue! 
19971 ) . Radial velocity surveys te nd to select stars 
that are chrom ospherically quiet ( Saar et al]|l998l : 
Gumming et a l. 2008), which is likely to introduce 
a bias toward older stars. Further, the link be- 
tween chromospheric activity and age means that 
distinguishing planetary transits from stellar sur- 
face features will pr obably be easi er for older, less 
active stars (but see lJenkin^ 2002 ). Consequently, 
there is a selection bias against planets around 
younger ("adolescent") stars. 

All of the solar system's "magnetic" planets 
(Earth, Jupiter, Saturn, Uranus, and Neptune) 
generate planetary-scale magnetic fields as the 
result of internal dynamo currents within the 
planet. The solar wind incident on these planetary 
magnetospheres is an energy source to the plan- 
etary magnetospheres, and the magnetosphere- 
solar wind interaction produces energetic (keV) 
electrons that then propagate along magnetic field 
lines into auroral regions, where electron cyclotron 
masers are produced. 

Specific details of the cyclotron maser emis- 
sion vary from planet to planet, depending upon 
such secondary effects as the planet's magnetic 
field topology. Nonetheless, applicable to all of 
the magnetic planets is a macroscopic relation 
relating the incident solar wind power Pgw, the 
planet's magnetic field strength, and the me- 
dian radio luminosity Lrad - Various investiga- 



Farrell et aLlll999t IZarka et alj l gOOll: iFarrell et al 

72 



dian radio luminosity l^rad- Various mvestiga- 
tors toesch fc Barrowlll984riDesch fc Kaiseilll984t 



Desch fc Ruckeijll985tlBarrow et al Ill986l : iRucker 
19871 : lDeschlll988l : iMillon fc Goertzlll98i ) find 



(1) 



with e the efficiency at which the solar wind power 
is converted to radio luminosity, and x « 1. The 
value for e depends on whether one considers the 
magnetic energy or kinetic energy, respectively, 
carried by the stellar wind. The strong solar wind 
dependence is reflected in the fact that the lumi- 
nosity of the Earth is larger than that of either 
Uranus or Neptune, even though their magnetic 
field are 10-50 times stronger than that of the 
Earth. 

The magnetospheric emissions from solar sys- 
tem planets and the discovery of extrasolar planets 
has motivated both theoretical ( Zarka et al.lll997 : 



2004':'Lazio et al.'2004':'StevenJ2005vGriessmeie r et al 
,2005; Zarka 2006, 2007; Gricssmcicr ct al. 2007a,^^ 



and observational work (iYantis et al.lll977l:IWinglee et al 



1986 : lBastian et al.|[2000l: Lazio et al.l2004lRvabov et al 



2004 : iGeorge fc StevensI l2007t iLazio fc Farrelll 
20071 : ISmith et al.ll2009h on magnetospheric emis- 
sions from extrasolar planets, including some be- 
fore the confirmed discovery of any extrasolar 
planets. 

Implicit in many of the early predictions for 
planetary magnetospheric emissions was that the 
stellar winds of other stars are comparable to the 
solar wind. Yet, from measurements of the sizes 
of astropauses (i.e., the boundary between the 
stellar wind and the loc al interstellar medium). 
Wood et all (|200l l2005f ) find the mass loss rate 
as a function of age, M oc t^, with x « —2, a de- 
pendence probably linked to the decrease in sur- 
face magnetic activity with stellar age. Thus, the 
stellar wind around a 1 Gyr old star may be 25 
times as intense as the current solar wind (from a 
4.5 Gyr old star). 

As a specific illustration of the possible effect 
of considerin g younger stars wi th more intense 
stellar winds, iFarrell et al. I (11999) and Lazio et al 



(|2004l ) predicted that the magnetospheric emission 
from the planet orbiting t Boo would be of or- 
der 1-3 mJy at frequencies around 30 MHz (10 m 
wavelength), assuming that the stellar wi nd of the 
star w as c omparable to that of the S un. IStevend 
(2005) and lGriessmeier et al. ( 2007a ). taking into 
account the likely stellar wind strength of t Boo, 
predicted that its emission would be at a level of 
order 100-300 mJy. The former prediction (1- 
3 mJy) is below the sensitivity of current instru- 
mentation; the latter is not. 

This paper reports a blind search for magneto- 
spheric emissions from nearby "adolescent" stars. 
In ^ we describe how we selected stars from ex- 
isting catalogs of nearby stars, in we present 
the results from our stacking analyses and as well 
as considering the stars individually, and in ij4]we 
present our conclusions. 

2. Stellar Catalog Assembly 

Three catalogs form the basis for our identifi- 
cation of nearby "adolescent" stars. 

1. The NStars program is part of the Space 
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Interfero metry Mission Prep arat ory Science 
Program. iGrav et all (|2n0.'^ and lGrav et aT 



( 20061 ) obtained spectroscopic observations 
of the 3600 main-sequence and giant stars 
within 40 pc of the Sun with spectral types 
ear her than MO. 



2. The Spectroscopic Properties of Cool Stars 
(SPOCS) is a compilation of stars form- 
ing the basis of radial velocity planetary 
searches. From the series of observations in 



the SPOCS program, iTakeda et al.l (|200<D 
estimated various physical parameters for 
the stars. 



3. The Geneva- Copenhagen survey of the solar 
neighborhood (CCS) is an effort to assem- 
ble a complete and consistent set of observa- 
tional and phy sical parameters for n earby F 
and G dwarfs ( Holmberg et al. lEooi). 



We applied four selection criteria — spectral 
type, distance, age, and declination — to the pub- 
lished catalogs to form samples for further study. 
For spectral type, we selected main-sequence F, 
G, or K dwarfs. The electron cyclotron maser — 
the same process by which planets generate ra- 
dio emission — has been detected from some stars. 
Typical detections are at frequencies of order 
5000 MHz, from which (lower) coronal magnetic 
field strengths of orde r 300 G or higher are inferred 



( Hallinan et al. l2008l) . For the electron maser 



emission process, the quantity B/v is approxi- 
mately constant, where B is the magnetic field in 
the emitting region. At frequencies ly « 100 MHz, 
we would expect magnetic field strengths of order 
30 G; lower frequencies would result in even lower 
field strengths. A field strength of 30 G is typical 
of that inferred for very late M dwarfs and early 
L brown dwarfs (Hcrgcr 2002) and within a fac- 
tor of a few of the field strength of Jupiter at the 
cloud tops (4 G). Moreover, while the Sun gener- 
ates intense emission at these frequencies, notably 
Type HI and IV radio bursts, even the strongest 
such radio bursts are far too faint to be detected 
over interstellar distances. Gergelvl (|l986) has 
considered the detection of solar-type stars at low 
radio frequencies and finds that current detec- 
tion thresholds are at least a factor of 10^ above 
what is required. If the targeted star is an F, G, 
or K dwarf, the presence of low radio frequency 



emission would be a strong indication that it is 
orbited by a sub-stellar companion. 

Two of the catalogs, NStars and SPOCS, con- 
tain only relatively nearby stars, with distances 
less than 40 pc and 130 pc, respectively; for the 
SPOCS catalog, fully two-thirds of the stars are 
within 40 pc. The GCS catalog extends to much 
larger distances, approximately 1000 pc. We im- 
pose a distance constaint for the following reason. 

In the general stacking analysis, as one stacks 
images from increasingly distant sources, the 
signal-to-noise ratio in a given image is decreasing, 
but is compensated by the increasing number of 
sources to stack. Considering stars to lie in shells 
at distances D with thickness AD, if all stars host 
planets, then the fiux densities of planets at larger 
distances will be smaller by a factor which 
will be exactly balanced by the increase in volume 
provided by going to larger distances, assuming 
a uniform distribution of stars. However, it is 
not yet known whether all stars host planets, and 
current limits are that the fraction of stars with 
planets is fp « 0.19, for plan ets with semi-majo r 
axes less than about 20 A\E jMarcv et al.ll2008h . 
Consequently, as one stacks images, an increasing 
fraction of the total number of images contain only 
noise so that the effective "signal-to-noise" ratio 
is decreasing. For the GSC catalog, we impose a 
distance constraint of 40 pc, for consistency with 
the other catalogs. 

Two of the catalogs (SPOCS and GCS) report 
the ages of the stars, which we adopt. For both 
catalogs, not only are ages reported but also an es- 
timated uncertainty. Given that stellar ages often 
have large uncertainties associated with them, in 
order not to exclude stars that might be younger 
than their nominal ages, we construct two sam- 
ples from both catalogs. One sample requires that 
the reported age of the star be less than 3 Gyr 
(samples that we denote SPOCS-age and GCS- 
age), while the other sample requires that the age 
taking its uncertainty account is less than 3 Gyr 
(samples that we denote SPOCS-eage and GCS- 
eage). The NStars catalog does not report the age 
of the star, but does report chromospheric flux 



' As magnetospheric emissions are powered by the stellar 
wind, for the purposes of this analysis, the presence of plan- 
ets more distant from their host star are not likely to be 
relevant. 



3 



in the Ca II H and K lines (logi?^^). We con- 
vert these to estim ated ages following the work of 



Henry et all (|l996h 



The final selection criterion, declination, is that 
we will focus on a northern hemisphere survey for 
which the effective declination limit is —25° ( ij3.ip . 
More generally, we are not aware of an all-sky cat- 
alog, at a common frequency and approximately 
constant rms image noise level, on which to per- 
form a similar analysis. Table [1] summarizes var- 
ious properties of the sub-samples from the three 
catalogs. 

This search is "blind" in the sense that we have 
made no effort to select stars already known to 
have planets, and even attempted to select stars 
whose properties are such that they might not pre- 
viously have been searched for planets. Nonethe- 
less, within our samples are small fractions of stars 
with known planets. We have cross-correlated our 
samples with the Extrasolar Planet Encyclopedia 
()Schneideil [2OO9I version July 1). The SPOCS- 
age sample has the highest fraction of stars with 
known planets, 6%, with the remaining samples 
all having fractions below 5%. The relatively 
large fraction of stars with known planets in the 
SPOCS-age sample is consistent with the larger 
SPOCS catalog being drawn from a set of stars 
which are or have been the focus of radial veloc- 
ity searches. The fraction of stars in our samples 
that are known to have planets is smaller than 
the estimated frac tion of stars with planets (19%, 
Marcv et al consistent with the notion that 

current planetary searches continue to be affected 
by selection effects. 

3. Magnetospheric Emission Searches 
3.1. Statistical Analysis 

Our search is b ased on the VLA L ow-frequency 
Survey (VLSS, ICohen et all 120071 ). This sur- 



vey imaged 95% the sky north of a declina- 
tion limit S > —30° at a frequency of 74 MHz 
(4 m wavelength) with a typical rms noise level 
of 100 mJybeam^^. Our focus on this survey 
is two-fold. First, the frequency of this survey 
is within a factor of 2 of the cutoff frequency of 
Jupiter (~ 40 MHz), so that it is plausible that 
other Jovian-like planets might emit at 74 MHz. 
Second, it is an electronically-available survey cov- 
ering a large fraction of the sky, so that many 



nearby stars are potentially accessible. There are 
other surveys at lower frequencies and which cover 
a significant fraction of the sky, to w hich this ap- 
proach might also be applied ( e.g.. iReedl 1990l : 

" '9m 



Dwarakanath &: Udava Shankail 119901 ) However 



the VLSS has the advantage of having images 
that are readily available in an electronic format 
combined with a high angular resolution and sen- 
sitivity. 

Each individual VLSS image was obtained by 
combining a series of "snapshots" acquired over 
a range of hour angles, with the time sampling 
within a snapshot being 10 seconds. A snapshot 
was typically 15-25 minutes in length, with snap- 
shots separated typically by about 1 hour. For 
comparison, Jovian decametric emission observed 
by the Nangay Decameter Array has been ob- 
served to have an average duration of about 1 hr, 
with a range from about 0.5 hr to a few hours, 
though these observations were probably domi- 
nated by the lo-contr oUed component of Jovian 



decametric emission ( Aubier et al. l2000l) . AL 



though it is too low in frequency to penetrate the 
Earth's ionosphere, Saturian kilometric radiation 
observed by the Cassini spa cecraft shows simi - 
lar temporal characteristics (jLamv et al. l2008h . 



but is not significantly affected by the presence 
of a major satellite. Assuming that extrasolar 
planetary magnetospheric emissions are similar to 
those of Jupiter and Saturn, if a planet was emit- 
ting during the course of a VLSS observation, it 
was likely to have been emitting for at least the 
duration of one snapshot, and potentially all of 
them. Consequently, there is a factor, potentially 
of order 1/3, in the luminosities that we derive be- 
low that would account for the fact that a planet 
might only be emitting for a fraction of the time 
that was used to acquire a VLSS observation. This 
factor is sufficiently small, relative to other uncer- 
tainties, that we shall not incorporate it explicitly 
into the analysis below. 

For each of our samples (Table [T]) , we down- 
loaded small images ("postage stamps", 1° in di- 
ameter) from the VLSS image server. Although 
the formal declination limit of the VLSS is —30°, 
the lowest declination fields are the most incom- 
plete and often have higher rms noise levels. Thus, 
we used an effective declination limit of —25°. 

We aligned each postage stamp image so that 
the target star was in the central pixel. The 
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Table 1 
Stellar Catalog Data 



Catalog 


Total Number 


Magnctosphcric 


Median 


Weighted Average 


Stacked Image 




Name 


in Catalog 


Subset 


Age 


Distance 


Noise Level 


Reference 








(Gyr) 


(pc) 


(mjy beam~^) 




NStars 


664+1676 


252+249 


1.3 


24.4 


5.7 


1,2 


SPOCS-age 


1074 


110 


1.9 


19.1 


11 


3 


SPOCS-eage 


1074 


176 


1.4 


14.3 


9.3 


3 


GCS-age 


16682 


355 


1.6 


21.4 


6.0 


4 


GCS-Eagc 


16682 


656 


0.7 


23.0 


4.8 


4 



Note. — The NStars catalog is published in two increments, a Northern Sample and a Southern Sample. Both the 
SPOCS and GCS catalogs provide age estimates and confidence intervals, from which we construct two measures of a 
star's age. The "age" samples adopt the nominal stellar age while the "cage" samples use the lower limit on the age 
estimate. Column 3, "magnetospheric subset," lists the number of stars from each catalog passing the four selection 
criteria of ^ 

References. — (DiGrav et al] ll2003t): (2')lGrav et al] ll2006h: (3')lTakeda et al.l ll2007t): (4) iHolmbere et al.l ||2009|) 



beam (point spread function) of the VLSS was 
80", with an image pixel size of 25". The coordi- 
nates of these stars are typically determined from 
the Hipparcos astrometric mission, epoch 1991.25 
(jPcrryman ct al. 1997); the VLSS observations 
were conducted between 2001 and 2007, with the 
majority of the observations conducte d between 
2003 September and 2005 April. Lazi o fc Farrell 
(j2007[ ) considered the possible astrometric uncer- 
tainties between the Hipparcos and VLSS frames, 
for the relatively high proper motion star r Boo 
(proper motion /i ~ 0"5 yr~^). They showed that 
the combination of astrometric uncertainties, un- 
corrected ionospheric refraction within the VLSS, 
and the proper motion of the star should have pro- 
duced an astrometric uncertainty in alignment of 
no more than 8", a fraction of a pixel. Thus, we 
are confident that alignment to the central pixel 
in the postage stamp images is sufficient. 

We examined each of the images for any sources 
that might be confused with a target star or that 
would be close enough to the location of a tar- 
get star to affect the stacking process. As an ex- 
ample, the star HD 69582, which appears in all 
of our samples, is approximately 75" (« 1 beam) 
from the radio source PKS 0814-029. This radio 
source is generally identified as a QSO, although 
no redshift has been measured. The offset between 
the star and radio source is significant enough that 



it is unlikely that the two are the same. We re- 
turn to the question of individual stars below. In 
addition, we found a small number of stars that 
are located close to the boundaries of the VLSS, 
particularly near the southern declination limit. 
The noise level in the images for these stars was 
much higher (and it is possible that the astrome- 
try is not as precise), so we excluded these. The 
total number of stars excluded on these consider- 
ations led to the sample sizes being about 10 stars 
smaller than a straightforward application of our 
initial selection criteria would suggest. 

The number of stars N in each of our sam- 
ples ranges from approximately 100 to several 
hundred. Assuming that the noise in the VLSS 
images is gaussian distributed, as generally ex- 
pected for radio intcrferometric imagesl^ we an- 
ticipate that the noise in a stacked image should 
be roughly 100A^~^/^ mJy beam"^, or about 5- 
10 mJy beam~^. In our stacking analysis, we com- 



^ In general, radio intcrferometric images of the sky are con- 
structed using a Fast Fourier Transform, with the flux den- 
sity of the visibility function at the spatial frequency origin 
taken to be identically zero. From Fourier Transform prop- 
erties, this so-called "zero-spacing" visibility value is equiv- 
alent to the total flux density within the field of view. In the 
absence of a source, the pixels in a thermal noise-limited 
image constructed in this manner will have a zero-mean 
normal distribution. The VLSS images that we analyze 
were constructed using this standard procedure. 
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bined the images in a weighted sense, weighting 
each image's contribution to the final stacked im- 
age by its individual rms noise level. Using the 
NStars sample as an illustration (Table [T|), and 
taking the rms noise levels in the images into ac- 
count, we expect that the rms noise level in the 
stacked image should be 5.6 mJy beam~^. The ac- 
tual rms noise level is 5.7 mJy beam~^, indicating 
that the assumption of gaussian noise-dominated 
images for the VLSS is justified, at least to the 
stacked image noise levels we have obtained here. 
Table [1] also presents the rms noise levels in the 
stacked images. 

In none of the stacked images do we detect sta- 
tistically significant emission. In an area of ap- 
proximately 1 beam in size, the strongest pixels 
range from approximately 1.5cr to 2.2cr. Table [2] 
presents 3a limits on the average fiux density of 
magnetospheric emissions from planets orbiting 
the stars in our samples from the various catalogs. 
Table [T] presents the weighted average distance for 
the stars in the various samples, with the weighted 
average distance for the N stars in a sample de- 
fined as 

4-e4^ (2) 



D2 



where Di is the distance to the i^^ star in the sam- 
ple. Using the weighted average distance, and as- 
suming that the bandwidth of planetary magneto- 
spheric emissions is comparable to the observation 
frequency, as it is in the case of Jupiter, we convert 
the fiux density upper limits to limits on average 
planetary luminosities. Table [2] also presents these 
planetary luminosity limits. 

3.2. Individual Stars 

Prior to stacking the images in each sample, 
we determined the peak intensity around each 
star relative to the individual image noise level 
(« lOOmJy beam"^). Stars for which the peak 
intensity exceeded 2.5tT were then re-examined. 

One motivation for performing this check is 
that the VLSS catalog of sources was constructed 
using a threshold test, relative to the image noise 
level, with a relatively high signal-to-noise thresh- 
old in order to maintain a low probability for 
a false detection. For the VLSS, the signal-to- 
noise thresholds was 7a. Thus, it is possible 
that there is stellar (or planetary) emission that 



would not have been cataloged. We found no stel- 
lar position for which radio emission above 3a 
(w 300 mJy beam~^) could be identified unam- 
biguously. There were stellar positions with radio 
emission above this level, but they could be ex- 
plained by other features, such as a sidelobe from 
another source. 

As noted above, the positions of some stars 
were close to, if not coincident with, radio sources, 
specifically the stars HD 38392, HD 49933, 
HD 79555, HD 143333, and HD 202575. These 



stars h ave be en detected by ROSAT (jHuensch et al 
1998 . 1999( ). with X-ray luminosities ranging 
from 6 X 10^" W to 3 x lO^^ W. The Benz-Giidel 
relation predicts that the centimeter-wavelength 
fiux densities of these stars should be of order 
0.1 mJy. Scaling these fiux densities to the VLSS 
(74 MHz, A = 4 m), we expect no emission to 
be detectable, for any reasonable radio spectral 
index. We also examined the NVSS (1400 MHz, 
20 cm) near the location of each star. The rms 
noise levels near these stars are comparable, and 
a 3a upper limit to the radio emission on any of 
these stars is 1.9 mJy, consistent with these stars 
not being radio sources, Finally, we are aware of 
targeted radio observations of only one of these 
stars, HD 143333, which placed the rather uncon - 
straining limit of 2 Jy at 5 GHz (jBlair et al.lll992| ) 
on the star or any associated planet. 



4. Discussion and Conclusions 

What do our limits imply about potential mag- 
netospheric emissions from planets orbiting any of 
the stars in our various samples? As noted in 
the discussion following equation the plan- 
etary luminosity can depend upon whether one 
is considering the kinetic energy Psw,kin or mag- 
netic energy Psw,mag carried by the solar wind. 
The kinetic energy power depends upon the stel- 
lar wind density n and velocity v as Psw,kin oc 
nv^ while the magnetic energy power depends 
upon the stellar wind velocity and magnetic field 



strength B as 



sw.ma.y 



cx vB^ ( Zarka et al. 2001 



Griessmeier et al. 2007b[ ). The strengths of all of 
these quantities both depend upon distance from 
the host star and are expected to evolve with stel- 
lar age. 



For t he specific case of t Boo, iGriessmeier et al 
(|2007bl ) illustrated how one can use a stellar wind 
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model for a star with a known age and a planet at 
a known orbital distance to estimate the strength 
of the planetary radio emission. A straightfor- 
ward extension of their approach could be ap- 
plied to stars not yet known to be orbited by a 
planet (s) and even samples such as those we have 
constructed here. The additional relevant quanti- 
ties that are needed would be the distribution of 
planetary semi-major axes and, in the case of an 
individual star, the distribution of its age estimate 
or, in the case of a sample of stars, the distribution 
of their age estimates. Combining these distribu- 
tions, one could estimate the appropriate scaling 
factor by which planetary magnetospheric emis- 
sion would be enhanced. However, in iJT]we argued 
that the current census of extragalactic planets is 
likely to be biased, particularly with respect to 
those planets that might be most likely to be ra- 
dio emitting. For that reason, we do not consider 
the distribution of planetary semi-major axes to 
be well enough constrained to incorporate it into 
our analysis, and we shall adopt a somewhat more 
simplified approach below. 

With respect to the stellar wind powers. 



Griessmeier et al.l (|2007al ) have synthesized var- 
ious observations and models to determine func- 
tional dependences for stellar wind velocity, den- 
sity, and magnetic field strength as a function of 
stellar age (their equations [15], [16], and [23]- 
[24]), for planets that are not too close to their 
host star. Applying these relations, and using the 
median age of the stars in the various samples (Ta- 
ble [T]), we determine a scaling factor, relative to 
the current solar value at 1 AU, for each of these 
quantities. From the scaling factors for the indi- 
vidual quantities (v, n, and B), the scaling factors 
for the kinetic energy and magnetic energy pow- 
ers are then determined. Table [2] (Columns 4-8) 
presents these scaling factors. 

These scaling factors are clearly only approxi- 
mate and somewhat model dependent. The num- 
ber of stars for which stellar wind parameters have 
been determined is small. Nonetheless, they serve 
as an indication of the potential effect of stellar 
age. We see that powers delivered (at 1 AU) to 
potential planetary magnetospheres around these 
stars, for most of these samples, may be enhanced 
by factors of 10-50 for kinetic power and by factors 
of 5-10 for magnetic power. 

The one potential exception is the GCS3-eage 



sample, for which much larger stellar wind amplifi- 
cation powers appear possible. These large factors 
result from the relatively small median age for this 
sample (0.7 Gyr, Tabled]). In turn, this small me- 
dian age likely reflects the relative lack of precision 
with which stellar ages can be determined. Many 
stars in the GCSS-eage sample have lower limits 
to their stellar ages around 0.1 Gyr, because their 
age estimates have large uncertainties (approach- 
ing 100%). 

These scaling factors for the stellar wind pow- 
ers are for a fiducial distance of 1 AU. As noted 
above, there is likely to be a strong star-planet 
distance dependence on the luminosity, but we do 
not attempt to include a distribution of planetary 
semi-major axes. Estimates of the star-planet dis- 
tance de2endencea£e_thatjt_js X < 2 
(e.g., iFarrell et all Il999t IZarkal [2003), implying 
that Jupiter would be about 25 times more lu- 
minuous were it at a distance of 1 AU instead of 
its current 5.2 AU distance. We therefore apply 
an additional scaling factor of 25. 

The final columns of Table [2] present Jupiter's 
luminosity scaled by these stellar wind kinetic en- 
ergy and magnetic energy factors, respectively, 
and the distance scaling factor of 25. Even if a 
planet with luminosity comparable to Jupiter or- 
bited every one of our target stars, our sensitiv- 
ity limits remain approximately a factor of 10-100 
above what would be needed to detect such plan- 
ets. If the typical planet-star distance is larger 
than our fiducial 1 AU value, the actual difference 
could be much larger. 

We have not attempted to combine the stacked 
images from the various samples. While the sam- 
ples themselves are homogeneous, they obtain age 
estimates from different methods. Further, there 
are many stars that are common to each sample, 
so that the stacked images are not independent. 

There are a number of next-generation, low 
radio frequency instruments under development. 
Notable among these are the Low Frequency Array 
(LOFAR) and the Long Wavelength Array (LWA). 
If they reach their design goals, both promise to 
provide rms noise levels a 3mJybeam^^, at 
frequencies below 100 MHz, nearly two orders of 
magnitude better than the 74-MHz VLA system 
which was used to conduct the VLSS. A simi- 
lar statistical analysis applied to future LOFAR or 
LWA observations may improve significantly upon 
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the limits presented here, or, ideally, detect extra- 
solar planetary emission. 

Should either LOFAR or the LWA detect emis- 
sion using a similar statistical analysis, identify- 
ing the stacked emission as planetary rather than 
stellar will be important. In targete d observations 



of an individual star or stars (e.g.. iBastian et al., 
2nnnl:lRvabov et al."2004V rxeorge fc Stevensll2nn7l : 



Lazio k Farrell 2QQ7; Smith et all l2009l ). an ob- 



vious distinguishing factor would be whether the 
emission is modulated with the planetary orbit. 
For this statistical analysis, alternate aspects of 
the stacked emission would have to be examined. 
Both LOFAR and the LWA are being designed 
to be broadband (over at least the 30-80 MHz 
frequency range), thus the radio spectrum of the 
stacked emission could be determined. Further, 
any correlation between the strength of stacked 
emission and the spectral type of the stars could 
be useful. 

Building upon LOFAR and the LWA will be 
the Square Kilometre Array (SKA) and the Lu- 
nar Radio Array (LRA). While their designs will 
be influenced by the work on LOFAR, LWA, and 
similar low radio frequency interferometers, both 
the SKA and LRA currently anticipate operating 
at frequencies that would be relevant for the de- 
tection of planetary magnetospheres; in the case of 
the SKA, the design goal for its lower operational 
frequency limit is 70 MHz, while, for the LRA, fre- 
quencies v ^ 50 MHz are envisioned. Both would 
likely provide an order of magnitude sensitivity 
improvement upon LOFAR and the LWA. 



We thank N. Kassim, S. Kulkarni, and B. Far- 
rell for discussions which spurred this analysis, 
A. Cohen for discussions on the VLSS, B. Erickson 
and B. Slee for discussions about solar decameter- 
wavelength emissions, and the referee for sev- 
eral comments that improved the presentation of 
these results. This research has made use of the 
SIMBAD database, operated at CDS, Strasbourg, 
France, and NASA's Astrophysics Data System. 
Basic research in radio astronomy at the NRL is 
supported by 6.1 Base funding. 
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Table 2 

Magnetospheric Emission Limits 



Sample 


Flux Density 
Upper Limit 
(mjy) 


Planetary 
Luminosity Limit^ 

(erg s~^) 


Velocity 


Density 


Stellar Wind Am 
Magnetic Field 


plification Factors 
Kinetic Energy 


Magnetic Energy 


K. E. Scaled Jovian 
Luminosity^ 

(org s"^) 


M. E. Scaled Jovian 
Luminosity^ 

(org s"^) 


NStars 


17 


9.0 X 10^^ 


1.7 


9.8 


2.4 


48 


9.5 


1.2 X 10^^ 


2.4 X 10^^ 


SPOCS-age 


33 


1.1 X 10^" 


1.4 


4.9 


1.8 


15 


4.8 


3.7 X 10^^ 


1.2 X 10^^ 


SPOCS-eage 


28 


5.1 X 10^^ 


1.6 


8.6 


2.2 


38 


8.3 


9.5 X 10^^ 


2.1 X 10^1 


GCS-age 


18 


7.3 X 10^^ 


1.6 


6.7 


2.0 


25 


6.5 


6.3 X 10^^ 


1.6 X 10^^ 


GCS-eage 


14 


6.8 X 10^^ 


2.2 


30 


3.6 


319 


28 


8.0 X 10^^ 


7.1 X 10^^ 



^1 W ^ 10^ erg s^^ 



Note. — The flux density upper limit and planetary luminosity limits are both 3cr. The flux density limits are converted to planetary luminosity limits using the weighted average distance 
for the samples (Table and assuming a 74 MHz bandwidth for the planetary magnetospheric emissions. We assume a fiducial planetary distance of 1 AU and a planetary luminosii 
scaling with distance as d^, so that the kinetic and magnetic energy scaled Jovian luminosities include an additional factor of 25; if the typical magnetospheric emitting planet is mo: 
distant from (closer) its host star than 1 AU, these Jovian scaled luminosities would have to be adjusted downward (upward). 
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